Influence of disorder on the perpendicular magnetoresistance of magnetic multilayers 
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The effect of disorder on the perpendicular magnetoresistance of magnetic multilayers is inves- 
tigated theoretically. Various kinds of disorder are considered: (i) interface substitutional disorder 
and (ii) bulk disorder in the various layers and in the leads. The calculations are based upon the 
non-equilibrium Green's function formalism, together with the recursion method for calculating the 
real-space Green's function. 
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I. INTRODUCTION 

The common wisdom about giant magnetoresistanoe 
(GMR) is that it is due to spin dependent scattering.!!] 
Although the transport is certainly diffusive for most sys- 
tems investigated experimentally, it has been pointed out 
that pure ballistic transport as well would lead to GMR 
for the perpendicular current geometry.0 Subsequently, 
it was proposed that for ballistic perpendicular trans- 
port, the electron confinement within the layers would 
yield a quantum size effect, i.e., an oscillatory behavior 
of the conductance and of the GMR as a function of layer 
thickness J§q 

In the present paper, we investigate the influence of 
disorder on perpendicular magneto-resistance, and dis- 
cuss the crossover from ballistic to diffusive transport as 
disorder increases. 



II. MODEL 

The model considered here consists of two ferromag- 
netic layers of thickness L separated by a non-magnetic 
spacer layer of thickness D; the structure is sandwiched 
between two ideal leads connected to ideal reservoirs. We 
use a simple-cubic single-band tight-binding Hamiltonian 
with a nearest-neighbor hopping. The hopping parame- 
ter, t, is the same for all pairs of nearest neighbors (for 
numerical calculations, t = — 1 will be assumed). The on- 
site energy of site i, £», depends on the chemical nature of 
the site and (for magnetic sites) on its spin. The majority 
spin on-site energy of the magnetic sites is taken equal to 
the on-site energy of nonmagnetic sites, and the minority 
spin on-site energy of the magnetic sites is 0.3 \t\ higher 
than the latter. The Fermi level is taken 0.8 \t\ above 
the bottom of the band; this corresponds to a quasi-free- 
electron case, with quasi-parabolic band structure and 
quasi-spherical Fermi surface. 

Two kinds of models for disorder are considered: (i) 
random substitutional disorder, and (ii) the Anderson 



model of disorder. Model (i) is particularly well suited 
to study the effect of interface interdiffusion. In model 
(ii), a random on-site energy (characterized by a square 
distribution of width W) is added to the on-site energy 
of the perfect case; this is a convenient model for bulk 
disorder. For both periodic (along the in-plane 

directions) model of disorder with 10 x 10 supercell is 
considered. We did not perform systematically an aver- 
age over disorder configurations, but we checked that this 
100-sites supercell is large enough for fluctuations to be 
unimportant. For k|| summations, we take 100 ky-points 
in the two-dimensional Brillouin zone corresponding to 
the supercell; this sampling is equivalent to taking 10 
k|| -points for the perfect case. 



III. THEORY 

The calculation of the conductance of the structure 
is based upon the non-equilibrium Green's function 
formalism^ When applied to the case of a two-terminal 
ballistic mesoscopic conductor we obtain the following 
result for the conductance:™] 



G = X Trk|l ' <T 



F; Gfi. T r G^) 



(3.1) 



In the above equation, the upper indices R and A re- 
fer to the retarded and advanced Green's function, re- 
spectively; the lower I and r indices refer to the left 
and right reservoirs, respectively. Thus, Gfjj: is the off- 
diagonal retarded Green's function linking the left reser- 
voir to the right reservoir; Ti (respectively, T r ) are given 
by IVr) — t 2 Au r ), where Ai (respectively, A r ) is the 
spectral density in the first plane of the left (respectively, 
right) reservoir, when it i s de coupled from the conductor. 
All the quantities in Eq. (3.1) are taken at the Fermi level. 
The conductance obtained here f rom the non-equilibrium 
Green's function formalism, Eq. (3.1), can be shown to be 



equivalent to the ones obtained from the Kubo formula, 
and from the Landaucr-Buttiker formalism.0 
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By separating the kii -conserving and k||-non- fusive contributions to the total conductance: 
conserving terms, we can separate the ballistic and dif- 



G ™ = T E Tr - [ T * ( k H ) G lr ( k ll ' k ll ) T r ( k ll ) G rl ( k ll . k ll )] . 
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(3.2a) 
(3.2b) 




disordered layer thickness (ML) 

FIG. 1. Conductance of a disordered nonmagnetic layer as 
a function of layer thickness. 



IV. RESULTS AND DISCUSSION 

Before considering the problem of perpendicular mag- 
netoresistance, we illustrate our method by presenting 
the conductance of a nonmagnetic layer of thickness up 
to 2000 monolayers (ML), with an Anderson-like disor- 
der of width W = 1.0|i|. The results are shown in Fig. [I]. 
For small thicknesses, the transport is quasi-ballistic, but 
for thicknesses larger than the mean-free path (approxi- 
mately 400 ML here), the transport becomes diffusive. 

Let us now consider the problem of the perpendicular, 
magneto-resistance. We define the magnetoresistance ascl 
A = (Gf — Gaf)/(Gf + Gaf)- The results are displayed 
in Fig. [| The magneto-resistance of the perfect system 
(Fig. ||a) exhibits the oscillations due to the quantum 
size effect .□ We have investigated the case of interface 
interdiffusion (Fig. ||b) by including two planes of ran- 
dom substitutional alloy: it appears that the latter influ- 
ences only weakly the magnetoresistance. We have also 
performed calculations treating the interface alloy within 
the virtual crystal approximation: the magnetoresistance 
(not shown here) is almost identical to the "exact" one 
(Fig. ||b); the effect of is interface interdiffusion is thus 
merely to replace the abrupt interface by a graded in- 



terface, which reduces the reflectivity, but produces little 
diffuse scattering. 

If we include bulk disorder in the spacer (Anderson 
model of disorder) the diffuse scattering leads to reduc- 
tion of the magneto-resistance for large spacer thickness, 
but the quantum size effect is still present (Fig. ||c). 

Next we simulate the effect of disoder in the leads by 
including a rather thick (400 ML) disordered region (An- 
derson model of disorder) in the left leado (Fig. gd) . This 
yields a significant reduction of magneto-resistance: this 
is easily explained as the effect of an additional resistance 
in series with the multilayer. The other effect is that the 
quantum size effect is strongly perturbed. 

In order to understand better the effect of disorder on 
the quantum size effect, we show in Fig. || the minority- 
spin contribution to the conductance of the ferromag- 
netic alignment (which is the one giving rise to the quan- 
tum size effect here). We see that, for a disordered lead 
(Fig ||d) , the amplitude of the conductance oscillation is 
reduced by a factor of 4, while the conductance itself is 
only reduced by 30 percent. Thus, this suggests that dif- 
fuse scattering in the leads supresses the quantum size 
effect (here it is not completely supressed because the 
disordered region is rather thin and the ballistic conduc- 
tance still accounts for 70 percent of the total conduc- 
tance). 
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FIG. 2. Magneto-resistance, (a) perfect case; (b) inter- Nagoya during which this work was initiated, 
face interdiffusion (2 layers of substitutional disorder); (c) 
bulk disorder in spacer layer (Anderson model, W = 0.5 \t\); 

(d) 400 disordered layers in left lead (Anderson model, 

W = 0.5 |*1 . 
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FIG. 3. Minority-spin conductance for the ferromagnetic 
configuration, (a)-(d): same as Fig. 
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